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1. 


La cinética heterogénea estudia la velocidad de las reacciones químicas en 


sistemas donde los reactantes se encuentran en diferentes fases o estados 


físicos, como sólidos, líquidos o gases. En estos sistemas, la velocidad de 


reacción depende de la interacción entre las moléculas de las distintas fases y 


la superficie de contacto entre ellas. 


Por otro lado, la cinética homogénea se refiere a las reacciones químicas que 


ocurren en sistemas donde los reactantes están en la misma fase, ya sea 


gaseosa, liquida o en solución. En este caso, las moléculas de los reactantes 


están mezcladas homogéneamente y pueden interactuar libremente entre si, 


facilitando el análisis de las velocidades de reacción. 


Aplicaciones de la cinética heterogénea en diferentes sistemas: 


Materiales: En el ámbito de los materiales, la cinética heterogénea es 
crucial para entender la corrosión de metales. La velocidad de corrosión 
de un metal esta influenciada por las reacciones electroquímicas en la 
interfaz entre el metal y el medio corrosivo. La formación de una capa de 
oxido en la superficie del metal puede afectar la velocidad de corrosión, 
actuando como una barrera. 

Procesos químicos: La cinética heterogénea resulta esencial en la 
catálisis. En las reacciones catalíticas, un catalizador solido interactúa 
con reactantes en fase gaseosa o liquida, acelerando la velocidad de 
reacción. La adsorción de los reactantes en la superficie del catalizador 
y la posterior formación de productor son etapas controlantes en este 
tipo de cinética. Ejemplos de aplicaciones incluyen la catálisis en la 
industria petroquímica, la síntesis de productos químicos y la 
purificación de gases de escape. 

Metalurgia extractiva: La cinética heterogénea es fundamental en los 
procesos de extracción y refinación de metales a partir de minerales. 
Durante la lixiviación, se utilizan reactivos químicos para disolver los 


minerales y extraer los metales. La velocidad de disolución depende de 
la superficie de los minerales y su interacción con los reactivos. La 
transferencia de masa, donde los iones metálicos difunden a través de 
una capa limite, puede ser una etapa limitante en este tipo de cinética. 


2. Simule la solubilidad de yeso y anhidrita en función de la temperatura. 
a) Reporte un gráfico [Ca+2] g/L - Temperatura (K) para cada uno. 

b) Gráfico dinámico densidad (kg/m3) - [Ca+2] g/L para tres temperaturas. 
Respuestas: 


e Primero se debe plantear el comando SOLUTION, para simular las 
propiedades de la solución acuosa, trabajaremos con agua pura, se 
indican los parámetros importantes como pH y temp, los otros se 
agregan por defecto. 

e Luego se introduce el comando EQUILIBRIUM_PHASE, el cual nos indica 
el equilibrio de fases y solubilidad de compuestos en solución acuosa, 
con Si = 0.0 que nos da condición en el punto de solubilidad o en 
equilibrio. 

e Luego el comando REACTION_TEMPERATURE el cual nos pide 
especificar el inicio el término de la temperatura para la reacción batch, 
este comando nos define el rango de temperatura en el cual ocurrirá la 
reacción. 

e Por último, el comando User. GRAPH para poder visualizar de manera 
grafica lo que ha sucedido en el proceso de disolución del Yeso (Gympse) 
y la Anhidrita (Anhydrite) 

e Hay que considerar que se utilizó la data base Phreeqc. 


Los resultados se y el código se muestran a continuación: 
Para el Yeso: 


SOLUTION 1 $ Agua Pura 
pH 7.0 charge 
temp 25.0 


EQUILIBRIUM PHASES 1 
Gypsum 0.0 10.0 


REACTION TEMPERATURE 1 
25.0 90.0 in 66 steps 


USER GRAPH 1 $ SOLUBILIDAD vs TEMPERATURA 

-chart title "Solubilidad en agua vs Temperatura, yeso" 
-headings Yeso 

-axis titles "Temperatura, (grados Celcius)" "[Ca], (g/L)" 
-axis scale x axis 25 90350 

-axis scale y axis auto auto auto 


-Start 
10 plot xy TC, tot("Ca")*40.08*TOT("water")/SOLN VOL, color = Blue, y axis = 1 
-end 
END 
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Para la Anhidrita: 


use SOLUTION 1 


EQUILIBRIUM_PHASES 1 
Anhydrite 0.0 10.0 


use REACTION TEMPERATURE 1 


USER GRAPH 1 # SOLUBILIDAD vs TEMPERATURA 

-chart title "Solubilidad en agua vs Temperatura, anhidrita" 

-headings anhidrita 

-axis titles "Temperatura, (grados Celcius)" "[Ca], ( g/L)" 

-axis scale x axis 25 90 5 0 

-axis scale y axis auto auto auto 

-start 

10 plot xy TC, tot("Ca")*40.08*TOT("water")/SOLN VOL, color = BLUE, y axis = 1 
-end 


END 


Solubilidad en agua vs Temperatura, anhidrita 
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Sin embargo, se pueden combinar estas dos solubilidades en un gráfico obteniéndose 


lo siguiente: 


SOLUTION 1 4 Agua Pura 
pH 7.0 charge 
temp 25.0 


EQUILIBRIUM_PHASES 1 


REACTION TEMPERATURE 1 
25.0 90.0 in 66 steps 


USER GRAPH 1 $ SOLUBILIDAD vs TEMPERATURA 

-chart title "Solubilidad en agua vs Temperatura, yeso" 

-headings Yeso 

-axis titles "Temperatura, (grados Celcius)" "[Ca], (g/L)" 

-axis scale x axis 25 90 5 0 

-axis scale y axis auto auto auto 

-start 

10 plot xy TC, τος ("Ca")*40.08*TOT ("water")/SOLN VOL, color = Blue, y axis = 1 
-end 


END 

use SOLUTION 1 

use EQUILIBRIUM PHASES 1 
use REACTION TEMPERATURE 1 


USER GRAPH 1 $ SOLUBILIDAD vs TEMPERATURA 

-chart title "Solubilidad en agua vs Temperatura, anhidrita" 

-headings anhidrita 

-axis titles "Temperatura, (grados Celcius)" "[Ca], (g/L)" 

-axis scale x axis 25 90 5 0 

-axis scale y axis auto auto auto 

-start 

10 plot xy TC, tot("Ca")*40.08*TOT("water")/SOLN VOL, color = BLUE, y axis = 1 
-end 


END 
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Respuesta (B): 


Para el gráfico dinámico se deben plantear 3 SOLUTION diferentes datos que ocurre a 
3temperaturas constante cero diferentes, escogimos 25, 50 y 75, el código se presenta 
a continuación el gráfico también: 


SOLUTION 1 $ Agua Pura 
pH 7.0 charge 
temp 25.0 


EQUILIBRIUM PHASES 1 
Gypsum 0.0 0.0 


REACTION 1 
Gypsum 1.0 
0.1 in 100 steps 


USER GRAPH 1 

-chart title "densidad del yeso, ϐ 25ος, 50ος y 75ος" 

-headings 25ος 50ος 75ος 

-axis titles "[Ca+2], (g/L)" "Density, (kg/m3)" 

-axis scale x axis auto auto auto 

-axis scale y axis auto auto auto 

-start 

10 plot xy tot("Ca")*40.08*TOT ("water")/SOLN VOL, RHO*(1/0.001), color = green 
-end 


END 


SOLUTION 1 # Agua Pura 
pH 7.0 charge 
temp 50.0 


use EQUILIBRIUM PHASES 1 
use REACTION 1 


USER GRAPH 1 

4 -chart title "densidad del yeso, 8 50ος" 

-axis titles "[Ca+2], (g/L)" "Density, (kg/m3)" 

-axis scale x axis auto auto auto 

-axis scale y axis auto auto auto 

-start 

10 plot xy tot("Ca")*40.08*TOT ("water")/SOLN VOL, RHO*(1/0.001), color - Blue, 
-end 


SOLUTION 1 £ Agua Pura 
pH 7.0 charge 
temp 75.0 


use EQUILIBRIUM PHASES 1 
use REACTION 1 


USER GRAPH 1 

* -chart title "densidad del yeso, @ 750C" 

-axis titles "[Ca+2], (g/L)" "Density, (kg/m3)" 

-axis scale x axis auto auto auto 

-axis scale y axis auto auto auto 

-start 

10 plot xy tot("Ca")*40.08*TOT("water")/SOLN VOL, RHO*(1/0.001), color - red, 
-end 


END 
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3. Simule la solubilidad de cal viva y cal hidratada en función de la temperatura. 
a) Reporte un gráfico [Ca+2] g/L - Temperatura (K) para cada uno. 
b) Gráfico dinámico conductividad (uS/cm) - [Ca+2] g/L para tres temperaturas. 


Respuestas: 


e Primero se debe plantear el comando SOLUTION, para simular las 
propiedades de la solución acuosa, trabajaremos con agua pura, se 
indican los parámetros importantes como pH y temperatura, los otros se 


agregan por defecto. 


Luego se introduce el comando EQUILIBRIUM_PHASE, el cual nos indica 
el equilibrio de fases y solubilidad de compuestos en solución acuosa, 
con Si = 0.0 que nos da condición en el punto de solubilidad o en 
equilibrio. 

Luego el comando REACTION_TEMPERATURE el cual nos pide 
especificar el inicio el termino de la temperatura para la reacción batch, 
este comando nos define el rango de temperatura en el cual ocurrirá la 
reacción. 

Por último, el comando User_GRAPH para poder visualizar de manera 
grafica lo que ha sucedido en el proceso de disolución de la cal viva 
(lime) y la cal hidratada (Portlandite) 

Hay que considerar que se utilizó la data base core.10 


Los resultados se y el código se muestran a continuación: 


Para Cal viva: 


SOLUTION 1 $ Agua Pura 
pH 7.0 charge 
temp 25.0 


EQUILIBRIUM PHASES 1 
Lime 0.0 10 


REACTION TEMPERATURE 1 
25.0 95.0 in 71 steps 


USER GRAPH 1 # SOLUBILIDAD vs TEMPERATURA cal viva 

-chart title "Solubilidad en agua vs Temperatura, Lime" 

-headings Lime portlandita 

-axis titles "Temperatura, (oC)" "[Ca], (g/L)" 

-axis scale x axis 25 95 5 0 

-axis scale y axis auto auto auto 

-start 

10 plot xy TC, TOT("Ca")*40.08/1000*TOT("water")/SOLN VOL, color = blue 
-end 


END 
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Para Cal Hidratada: 


use SOLUTION 1 


EQUILIBRIUM_PHASES 2 
Portlandite 0.0 10 


use REACTION TEMPERATURE 1 


USER GRAPH 2 $ SOLUBILIDAD vs TEMPERATURA Cal hidratada 

-chart title "Solubilidad en agua vs Temperatura, Cal Hidratada" 
-headings Portlandite 

-axis titles "Temperatura, (grados Celcius)" "[Ca], (g/L)" 

-axis scale x axis 25 95 5 0 

-axis scale y axis auto auto auto 

-start 

10 plot xy TC, TOT ("Ca") *40.08*TOT ("water")/SOLN VOL, color - red 
-end 


END 


Solubilidad en agua vs Temperatura, Cal Hidratada 
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Ahora bien, también como se realizó anteriormente se puede visualizar el 
comportamiento de estas dos especies en el siguiente gráfico: 


SOLUTION 1 $ Agua Pura 
pH 7.0 charge 
temp 25.0 


EQUILIBRIUM PHASES 1 
Lime 0.0 10 


REACTION TEMPERATURE 1 
25.0 95.0 in 71 steps 


USER GRAPH 1 # SOLUBILIDAD vs TEMPERATURA cal viva 

-chart title "Solubilidad en agua vs Temperatura, Lime" 

-headings Lime portlandita 

-axis titles "Temperatura, (oC)" "[Ca], (g/L)" 

-axis scale x axis 25 95 5 0 

-axis scale y axis auto auto auto 

-start 

10 plot xy TC, TOT("Ca")*40.08/1000*TOT("water")/SOLN VOL, color = blue 
-end 


END 
use SOLUTION 1 


EQUILIBRIUM PHASES 2 
Portlandite 0.0 10 


use REACTION TEMPERATURE 1 


USER GRAPH 1 $ SOLUBILIDAD vs TEMPERATURA Cal hidratada 

$ -chart title "Solubilidad en agua vs Temperatura, Cal Hidratada" 
$ -headings Portlandite 

-axis titles "Temperatura, (grados Celcius)" "[Ca], (g/L)" 

-axis scale x axis 25 9550 

-axis scale y axis auto auto auto 

-start 

10 plot xy TC, TOT("Ca")*40.08*TOT("water")/SOLN VOL, color = red 
-end 


END 
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Respuesta (B): 


Para la conductancia se logro crear el siguiente código: 


SOLUTION 1 H20 pura 
temp 20 
pH 7 


REACTION 
Portlandite 10 
5 in 100 steps 


EQUILIBRIUM_PHASES 
Portlandite 0.0 0.0 


USER_GRAPH 1 

-chart title "propiedades de Ca-H20" 

-headings 250C 30ος 40ος 

-axis titles "Ca g/L" "Specific conductance " 

-axis scale x axis auto auto auto 

-axis scale y axis auto auto auto 

-start 

10 plot xy tot("Ca")*40.08*TOT("water")/SOLN VOL, SC, color = green, y axis = 1 
-end 


END 


SOLUTION 2 $ Agua Pura 
pH 7.0 charge 
temp 30 


use REACTION 1 
use EQUILIBRIUM_PHASES 1 


USER_GRAPH 1 

+ -chart title "Conductividad de la cal hidratada, ϐ 30ος" 

-axis_titles "Ca g/L" "Specific conductance " 

-axis scale x axis auto auto auto 

-axis scale y axis auto auto auto 

-start 

10 plot xy tot("Ca")*40.08*TOT ("water")/SOLN VOL, SC, color = red, y axis = 1 
-end 


END 


SOLUTION 3 # Agua Pura 
pH 7.0 charge 
temp 40 


use REACTION 1 
use EQUILIBRIUM PHASES 1 


USER GRAPH 1 

+ -chart title "Conductividad de la cal hidratada, @ 40ος" 

-axis titles "Ca g/L" "Specific conductance " 

-axis scale x axis auto auto auto 

-axis scale y axis auto auto auto 

-start 

10 plot xy tot("Ca")*40.08*TOT ("water")/SOLN VOL, SC, color = Blue, y axis = 1 
-end 


END 


Donde se obtuvo el siguiente resultado: 


propiedades de Ca-H20 
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Specific conductance 


Sin embargo, esto nos causa duda, dado que a pesar de que la cal sea muy poco 
conductora, por los ¡ones disuelto de la solución debiera existir al menos un poco de 
conductividad en el gráfico, sin embargo, esto no es así, el posible error es que haya 
fallado el código, sin embargo, se probó el mismo código para la halita en las mismas 
condiciones de temperatura y este si arrojó resultados correctos para tal código, se 
muestra a continuación el resultado del mismo código presentado: 
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8. (G5, G6) Simule la cinética de disolución de Pirita en Phreeqc. Realice un análisis de 
sensibilidad al graficar la concentración de [Fe] mg/L -- tiempo - pH. Comente sus 
resultados. 


a. Disolución en agua pura, ARD y PLS 


SOLUTION 1 8 disolución en ARD SOLUTION 1 + PLS 
SOLUTION 1 $ agua Pura SIRITEFON CARD units mol/L 
Ν temp 25 
units mol/L units mg/L temp 25 
temp 25 redox pe pH 1:3 
7 pH 2.2 
pH dd Cu+2 0.38 
K 50 S04-2 0.9 
Ca 230 
ls 1200 Cl- 0.98 
Cl 550 charge Fe+2 7.16e-3 
C(4) 4 
S(6) 13500 Fe+3 0.02 
Fe (3) 1500 Μπ΄ «2 0.03 
Al 800 
A1+3 0.2 
Mg+2 0.4 
AAA AAA } PYRITE ahopra se debe simular la pirita con los datos RATES DE LA PÁGINA 
RATES $ este es el modelo cinético para la pirita 
Pyrite dissolution $22 ponerle nombres a las variables 
-start 
rem PARM(1) = Time (s) 
rem PARM(2) = Temperature (K) 
10 Ea = 5€000 $J/mol 
20 Ec = 0 $J/mol 
30 ta = 1.578e7 $$ segundos 
40 MD = 7.87 $ densidad del hierro 


50 Atwt = 55.845 $ peso molecular Fe 
€0 ka = 1.58548e-12 $ cm/s 

70 Tref = 298 $ en kelvin 

80 Kc = 1.585485e-14 $ cm/s 

9508 = 8.3144€21 $ ctante de los gases 
100 SA = 50300 $área de la superficie 
110 if (m <= 0) then goto 1€0 

120 ratel = -Ka*EXP(-PARM(1)/ta)*EXP[(-Ea/R) *(1/PARM(2) - l/Tref)]-Kc*l*EXP[(-Ec/R)*(1/PARM(2) -1/Tref)] $ en cm /s 
130 rate2 = ratel*(MD/Atwt) $mol/cm2/s 

140 rate = rate2*SA $mol/s 

150 moles = rate*TIME 

1€0 SAVE moles 


KINETICS 8 una vez completado el modelo cinético se debe utilizar este comando para definirlo 
Pyrite dissolution 


-m 0.14 $mol inicial 

-formula Fe -1 $ se pone la fórmula del elemento 

-parms 8.888e€ 298 $ parm l y parm 2 para que ingrese en el RAtes 
-step 8.888e€ in 100 steps 

-step_divide 100 pasos 

-runge_kutta € $ número de integración 

-tol le-15 $ tolerancia error 


INCREMENTAL REACTIONS true 


USER_GRAPH 1 

-chart title "Análisis de sensibilidad" 

-headings time Fe pH 

-axis titles "Tiempo, (días)" "Fe, (mg/L)" "pH" 
10 GRAPH x TOTAL TIME/8€400 

20 GRAPH y TOT("Fe")*55.85e3 4 multiplicación lo deja en mg/L 
30 GRAPH SY -LA("H+") $ esto es el pH 


-end 


END 


Para este ejercicio se planteó el RATES y KINETICS correspondiente y luego se 
cambiaron las soluciones para notar el efecto en las gráficas. 
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Análisis de sensibilidad en PLS 
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Se puede observar que el comportamiento del pH varía dependiendo de las distintas 
soluciones en las que esté expuesto, para la solución de agua pura se observa que el 
pH tiende a incrementar pero de manera no tan abrupta, en cambio para la solución en 
PLS se puede observar que se forma un escalón con un crecimiento del pH alto en poco 
tiempo entre los 20 y 40 días, por otro lado está la solución ARD y este se comporta de 
manera escalonada teniendo cambios fuertes de pH entre los días 20,40 y en los días 
75 y 80. 


Se pude observar que entre menos cambios escalonados tenga el pH en la solución 
más alto será su valor de pH, ya que la solución de agua pura no tuve un 
comportamiento escalonado y se obtuvo un pH de aproximadamente 13. 


b. Evaluar cinética a tres temperaturas a elección 


Para este caso se utilizaron las temperaturas 298,323 y 473 para cada solución. 


Solución Agua pura 
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Análisis de sensibilidad a 473 kelvin 
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Se puede observar cómo la temperatura afecta el incremento del pH, si bien a mayor 
temperatura aumenta el pH en menores tiempos, hay que considerar también que el 
máximo es menor a medida que se aumenta la temperatura. 


El máximo disminuyó cada vez que se aumentaba la temperatura, a la temperatura de 
298 kelvin el pH se acercaba a 13, en cambio cuando la temperatura era 473 el pH bajó 
aproximadamente hasta 8,4. 


c. Diferentes presiones parciales de oxígeno en solución. 


Para este caso se agregaron diferentes concentraciones de O a la solución agua pura 
para ver el efecto que tiene. 
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11. (G5, G6) Simule la cinética de disolución de Quarzo en Phreeqc. Realice un análisis 
de sensibilidad al graficar la concentración de [Si] mg/L - tiempo. Comente sus 
resultados. 


a. Disolución en ADR, PLS y agua de mar 


SOLUTION 1 48 Agua de Mar SOLUTION 1 $ PLS SOLUTION 1 ARD 
units g/L units mol/L temp 25 
temp 25 $ oC " 
f temp 25 units mg/L 
> — pH 1.3 redox pe 
Na 10.579 Cu+2 0.38 pH 2.2 
Mg 1.21 504-2 0.9 Na 400 
Ca 0.4102 charge Cl- 0.98 K 50 
κ 0.386 Fe+2 7.16e-3 Ca 230 
cl 19.492 Fe+3 0.02 Mg 1200 
S(6) 2.552 as SO4 Mn'«2 0.03 61 550 charge 
C(4) 0.119 as HCO3 A143 0.2 C(4) E 
Mg+2 0.4 S(6) 13500 
Fe (3) 1500 
Al 800 


Se agregaron las diferentes soluciones al mismo código sin variar RATES Y KINETICS. 


RATES £$Modelo cinético 

Quartz fEspecie de interes 

-start 

10 moles-(parm(1)/parm(2))*((m/m0)^(2/3))*(10^-13.7)* (1-SR("Quartz")) 
20 save moles*time 

-end 


KINETICS Parámetros del modelo cinético 

Quartz 

-formula 5102 

-m0 102.7 £moles iniciales 

-parms 22.7 0.162 £$Area superficial,Porosidad 

-step 1.58e8 in 10 steps £Segundos 

-tol le-8 fTolerancia intervalo de integración Solver Runge Kutta. 
INCREMENTAL REACTIONS true 


USER_GRAPH 1 

-chart_title "Cinética de disolución quarzo" 
-headings Tiempo Si 

-axis title "Tiempo (Años)" "Si (mg/L) " 
-start 

lO graph x Total time/3.1536e7 

20 graph y TOT ("Si")*28.08e3 

-end 

END 


Cinética de disolución quarzo en agua de mar 
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Cinética de disolución quarzo en PLS 
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Cinética de disolución quarzo en ARD 
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Se puede observar que para los tres casos el 


máximo de pH son distintos y que la 
solución ARD destaca por tener una 
cinética más rápida que las otras 
soluciones, pero la solución que logra tener 
un mayor pH es la de agua de mar, esto se 
presenta por los diferentes elementos 
químicos que existen en las soluciones y 
cómo afectan en la cinética. 


b. Evaluar cinética a tres temperaturas a elección. 


Para este caso se utilizó la solución agua de mar y se analizó como afectaba la 
cinética a 100,200 y 2? C 


Cinética de disolución quarzo en agua de mar a 100 *C 
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Cinética de disolución quarzo en agua de mar a 2 °C 
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Cinética de disolución quarzo en agua de mar a 200 *C 
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Se puede observar como la 
cinética a mayores 
temperaturas llega a un valor 
más alto y en menos tiempo, 
en cambio si la temperatura es 
baja la cinética es más lenta y 


llega a concentraciones más 


bajas. 


